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COMPUTER CONTROLLED ANTENNA SYSTEM 
Nicholas A. Raumann 


ABSTRACT 

With the advent of relatively cheap and small computers on the 
market, digital techniques become a very attractive application to the 
servo and control system of large antennas. A small dedicated digital 
computer can be used economically to replace certain pieces of existing 
hardware and can be used, in addition, to perform an automatic checkout 
and readiness verification of the servo system. Such a system has several 
advantages: reliability and consistency of operation, minimization of turn 
around time between satellite passes and fast identification of component 
or subsystem failures or misadjustments. 

These techniques have been evaluated at the Network Test and Train- 
ing Facility (NTTF) using the 40-foot antenna and the Sigma V computer. 
Programs have been completed which drive the antenna directly without 
the need for a servo amplifier, antenna position programmer or a scan 
generator. In addition, several programs have been developed which test 
the operation of the servo system and diagnose subsystem failure or 
misadjustments. 
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COMPUTER CONTROLLED ANTENNA SYSTEM 


1. INTRODUCTION 

With the advent of relatively cheap and small computers on the market, 
digital techniques become a very attractive application to the servo and control 
system of large antennas. A small dedicated computer can be used to perform 
automatic checkout and readiness verification of the servo system. In addition, 
the computer can be used to replace certain pieces of existing hardware. The 
tracking performance of the antenna can also be improved by use of improved 
control algorithms. 

| These techniques have been evaluated at Goddard’s Network Test and Training 
Facility using the 40-foot data acquisition antenna and an XDS/Sigma V computer. 
Even though this computer is a much larger machine than is necessary for this 
application, the feasibility of the various techniques has been established. All 
real time programs have been written in assembly language (SYMBOL), while 
most test programs have been written in FORTRAN IV for reasons of expediency. 
Programs have been evaluated under realistic conditions and system operation 
has been demonstrated by tracking NIMBUS and OGO satellites. 

Software developed in support of the monitoring and automatic checkout 
function is shown in Figure 1.1. During normal operation of the antenna system 
the computer monitors critical quantities such as power supplies, oil tempera- 
tures, axis velocities and limits. Abnormal conditions are flagged to the 
operator and if condition warrants it, direct action is taken by the computer. 

The antenna operator can call for an initialization program, which automatically 
brings the antenna to operational status by properly sequencing the various 
switching functions. 

A readiness verification program can be called which is intended as a short 
prepass check of the servo system. This program checks power supplies, 
circuit breakers, interlocks, and the hydraulic system. As a final test a step 
function is applied to the antenna and its transient response is evaluated according 
to a predetermined performance index. If all of these tests pass satisfactorily, 
the antenna is declared operational. If not, the operator is informed of the 
particular problem and is given the choice of operating the antenna in a degraded 
mode or calling the complete diagnostic test program which will help in further 
pinpointing the problem. 
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Figure 1-1. 


Programs for Monitoring and Automatic Checkout of Antenna. 


The diagnostic test program checks for items such as tachometer gain and 
ripple, error sensitivity, synchro operation, friction, low speed performance, 
and hydraulic leakage; and finally generates a diagnostic message to the operator. 
This program can also be called during the performance of routine maintenance, 
because many tests of this program coincide with tests required by standard 
maintenance procedures. 

Figure 1-2 shows the tracking loop for an antenna as it exists at a typical 
STDN tracking site. The antenna structure is driven by a hydraulic drive sys- 
tem, while the tachometer, receiver and encoder provide feedback signals from 
the structure. The tracking loop is closed via the servo amplifier and the 
receiver or the antenna position programmer depending on the operating mode 
selected. Off-line the PB 250 computer (UNIVAC 1218 for USB antennas) is used 
to generate predict data and a scan generator is provided which aids during the 
satellite acquisition phase. 

All devices shown crosshatched in Figure 1.2 are special purpose equipment 
ranging in size from one chassis to several racks. The function of these devices 
have been replaced with a program on the same computer. Test results indicate 
that the digital system has a performance that is equal to or better than the 
performance of the existing system. 
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Figure 1-2. Tracking Loop of Antenna. 


Using the computer approach, the tracking performance of the antenna can 
be improved. This has been accomplished by addition of a new mode of operation, 
termed "augmentation mode." This mode permits operation in the program 
mode with the addition of a correction factor derived from the receiver and 
passing it through a low-pass filter and integrator. Thus, the antenna is driven 
by a signal with high signal- to- noise ratio, as it would normally be during pro- 
gram mode operation but improved pointing capability is provided by the intro- 
duction of the average receiver output. 

Conceptually, this project has been outlined in reference 1.1. 
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2. MONITORING AND AUTOMATIC CHECKOUT 


The Monitor, Initialization, Readiness Verification and Diagnostic Test 
programs identified in Figure 1.1 consist of a calling sequence of various sub- 
routines which test or measure certain aspects of the antenna aystem. These 
programs form the link between the computer, the antenna system and the 
operator and since they consist of logical manipulations of subroutines they will 
not be discussed in detail. The discussions below will concentrate on techniques 
employed in the various subroutines and is largely based on a study covered in 
reference 2.1. In these discussions reference is made to various operating 
modes of the antenna and figure 2.1 is included to show the implementation of 
these modes in the antenna system. 

2.1 AC Power 

In this application various amplitudes have to be measured at a nominal 
frequency of 60 Hz. Amplitude is determined by taking 300 samples at a 12 KHz 
rate. From these samples, 200 consecutive samples are chosen, such that the 
first sample has a magnitude of 10% or less of the absolute peak value. If such 
a sample complement cannot be found, all 300 samples are rejected, and a new 
sampling process is initiated. This procedure insures that the measurement 
accuracy will be below a certain value and at the same time the sampling process 
will take less than 25 ms. 

The root -me an- square (rms) value of the voltage rms is then determined 
by using: 


V 


rms 


i 

200 


200 



1/2 


E 


2 

K 


where E k = voltage sampled at k-th interval, volts 

The accuracy of this measurement including the analog-to-digital converter 
accuracy is ±0.1%. 

The frequency is established by scanning the 300 samples and determining 
the number of samples, N, from polarity change tc polarity change of the same 
sense. Frequency f in Hz is then computed using 

f _ 12000 
N 

The accuracy of this measurement is ±0.5%. 
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Figure 2-1. Simplified Diagram of the Control and Drive System 













2.2 DC Power and Hydraulic Pressure 


DC voltage levels are measured by taking 200 samples at 12 KHz rate. The 
average value of the voltage, V ave , is computed by 


V 


ave 


200 



K=1 


( 2 . 1 ) 


This measurement is accurate to ±0.1%. 

Ripple, V rms , riding on this dc voltage is established using: 


V 

rms 



( 2 . 2 ) 


This measurement is accurate to ±2%. 

Hydraulic pressure is measured by means of pressure transducers whose 
output is a dc voltage proportional to pressure. Hence, measurement involves 
determination of V ave and V fms using equations (2.1) and (2.2), and applying 
an appropriate conversion factor. 

2.3 Tachometer Gain and Ripple 

In this test a constant voltage is applied to the rate loop of the antenna. 
After the decay of transients, tachometer output and encoder are sampled at a 
100 Hz rate until 500 sample pairs have been collected. The average and rms 
tachometer voltages are computed using relationships (2.1) and (2.2) only 
changing the summation limit to 500. Antenna velocity, q , is established 
from the encoder readings: 


e 

ave 


i 

1250 


500 




(2.3) 


where 


6 - encoder reading at k’th interval. 

K 

& Q = initial encoder reading. 
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From this tachometer gain is computed: 


Gain 


ave 

e 


ave 


Ripple is computed by establishing first the velocity ripple, <9 rms , of the antenna 
as measured by the encoder and then forming the difference in root-mean- square 
fashion between rms tachometer voltage and 6 : 

0 rms 


6 = 
rms 



Ripple 


V \ 

rms l 

.gainJ 


- 6 2 

rms 


1/2 


(2.4) 


2.4 Low Rate Performance 

This test is used to establish whether the antenna can operate smoothly at 
a low velocity and is performed for the rate loop and then repeated for the 
position loop. A constant velocity input is provided for either loop and after 
decay of transients, 500 encoder samples are taken at a 100 Hz rate. Velocity 
and ripple are computed using relationships (2.3) and (2.4). 


2.5 Deceleration Test 

This test is performed to assure that maximum antenna deceleration is 
within certain limits and that the antenna comes to a complete stop within a 
specified time. The antenna is moved at a constant velocity and the brakes are 
applied. Prior to application of brakes 500 encoder and 500 tachometer samples 
are taken to a 100 Hz rate. From this the average velocity 6 ave and the 
average tachometer voltage V ave are computed using equations 2.3 and 2.1 
respectively. After the application of the brakes 1000 tachometer voltage 
samples are taken at a 800 Hz rate. Deceleration, 9 , is computed by means 
of the following expression: 
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1600 


e = 


(K b -K y V av 


6 v — i 

ave \ 


(V k - V a ) 


K=Ka 


where K a = 1st sample at which V k < 0.9 V 

= 1st sample at which V ^ 0.1 V 


After establishment of 8 , the remaining tachometer voltage samples are scanned 
and the total stopping time is determined. 


2.6 Breakaway Torque 


This test measures the breakaway torque and hence, the amount of static 
friction encountered in the antenna system. A step function of sufficient magni- 
tude is applied to the rate loop., 500 samples pairs are taken of the differential 
pressure, P, across the hydraulic motors and of the encoder at a rate of 100 Hz. 
The sample, K b , is determined at which the encoder changes value. Thus the 
breakaway torque, T, is computed: 


T = 



where K m = motor torque constant. 


2.7 Error Sensitivity 


This test is used to establish the error sensitivities in the various positional 
modes of operation of the antenna. A particular mode is selected and the antenna 
is permitted to null out. Then 200 initial sample pairs are taken of the error 
voltage, e , and encoder output at a rate of 100 Hz. The initial averages are 
computed for the error voltage, e i ave , and encoder, d i ave , using equation (2.1). 
Then an offset voltage of approximately 1 degree is introduced. After transients 
have died away another 200 sample pairs are taken and final values for error 
voltage, e f ave , and encoder d f avg are computed. The error sensitivity, S, 
is then established: 
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(2.5) 


S = 


6 


f ave 


€. 

1 ave 


e c - e. 

f ave i 


ave 


The amount of noise on the error signal is then calculated using equation (2.2) 
to find e. rms and e f rms . The total noise is then defined as: 


Noise 


1 [e? 

2 i rms 


e2 f ] 

f rms J 


1/2 


The above procedure is repeated for all positional modes. 


2.8 Synchro Alignment 

Synchros are used in the manual position, slave and stow modes. The 
alignment of these synchros is checked by entering the various modes and read- 
ing the encoders and error voltages after transients have disappeared. 200 
sample pairs are taken at a 100 Hz rate and the average encoder reading, d , 
and the average error voltage, e gve , are determined using equation (2.1). The 
synchro alignment, 6 S , is then found: 


e 


= e 

ave 


+ S e 


ave 


where the error sensitivity, S, has been established in equation (2.5). 


2.9 Hydraulic Drive System 


The schematic of the hydraulic drive system is shown in Figure 2.2 (Only 
one axis drive is shown.) In this test the antenna is operated at a constant speed 
of 4°/sec and then l°/sec. For both speeds pressures P x , P 2 , P 3 , P^ and p ; 
differential pressures A P t and A P 2 ; flow rates Q } and Q ; and valve current I v 
are measured. Measurement follows the procedure discussed in Section 2.2. 

First, the filter constants K fl and K f2 are computed which represent the 
degree to which the filters are clogged (using 4°/sec data): 


K fi = 


A^ 

Qi 


and K f 2 = • 


AP„ 


Q 2 
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Then the line drops K el and K e2 are established (using 4°/sec data): 



and 



Subsequently, drive system leakage is computed by using the following equation 


Qi = K p 3 + L d (P 3 


- P 4 ) + D 


m 


e 

m 


where 


L m = internal motor leakage 

L, = leakage through damping valve 

d 

D m = displacement constant of hydraulic motor 
6 = motor shaft velocity 

m 

By substituting values into above equation corresponding to antenna speed of 
4°/sec and then l°/sec, the two unknowns L m and L d can be numerically 
evaluated. 

Finally, the gain of the servo valve, K s , and the gain of the hydraulic drive 
system, K h , are established: 


I v/p - P 
i v r 2 3 

and 

6 

m 


Equations in this section have assumed that motor number one is driving or 
P 3 > P 4 . Equivalent analysis can be performed on motor number two. 

2.10 Bate Loop Performance Index 

The rate loop performance index is a measure of the transient response of 
the rate loop; a lower number of the index indicating a better performance. The 
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test is performed by applying a step function to the rate loop and sampling the 
velocity error, e v , at a rate of 100 Hz until 500 samples have been collected. 
The rate loop performance index, RPI, is computed using a formula that was 
empirically developed: 


RPI 


RT 

0.07 


+ 


10 


e 

V rms 

> 

e 

v ave 


RT :> 0.07 


where 


RT = rise time, s 

e v rms ’ e v ave ~ stea( *y state ripple and average velocity error respectively 
using last 300 samples. 

The characteristics of such a performance index can be judged from Figure 2.3 
where the index is plotted vs. loop gain. Also shown in that figure are step 
function responses for low, nominal and high loop gain. 


2.11 Position Loop Performance Index 

This performance index is essentially similar to the previous one except 
that the transient response of the position loop is evaluated. 1000 samples of 
the position error are taken at a 100 Hz rate. These samples are taken for both 
type I and type II operation of the position loop. (Type n operation corresponds 
to a loop having an additional integration, whereas type I does not have this 
integration.) Empirically the position loop performance index for type I opera- 
tion, PPIj, has been established as: 


PPI = J_ (2^ + ST); ST £0.75 
1 0.75 \ 8 / 

where 

OS = overshoot, percent 
ST= settling time, seconds 

(= elapsed time up to a point after which response stays within 
5% of final value.) 

The index for type II operation, PPI 2> is: 
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OS :> 40, ST £ 1.5 


PPI 2 = — 
2 3.5 



ST\ 
3 /’ 


Graphically the characteristics of these indices are shown in Figures 2.4 and 2.5. 

It should be noted here that experimentation has been done with more con- 
ventional indices such as ISE (integral of the square error), ITAE (integral of 
product of time and absolute error) and several others. However, all of these 
indices showed a lack of sensitivity around the nominal setting. 


3. TRACKING LOOP 

As was pointed out in the introduction the tracking loop can be closed via the 
computer replacing several pieces of special, purpose hardware. In this project 
an attempt has been made to replace hardware functions on a one-to-one basis, i.e. 
no improvements or optimizations of the tracking loop have been attempted. 
However, it should be stressed, that if improvements or changes are required, 
the computer approach yields a high degree of flexibility. 

The antenna system, as it exists in the field today, is operated from a control 
console which is an arrangement of several racks. The operator is provided 
with several displays and controls which permit complete operation of the antenna. 
A low cost computer driven CRT display with an associated keyboard and some 
special control functions can be used to replace the present control console. It 
is anticipated that such an approach will enhance human engineering aspects of 
antenna control and will also represent considerable savings when antenna pro- 
curements are considered. The CRT display approach is under investigation, 
but is not reported here. 

Details of the program structure and techniques are covered in reference 3.1. 

3.1 Digital Compensation 

The tracking loop in all modes of operation is closed via the servo amplifier. 
Inputs are command signals and signals from feedback transducers while the 
output is a voltage which drives the valve of the hydraulic drive system. In the 
servo amplifier the drive signal is shaped in accordance with stability and track- 
ing accuracy requirements. This signal shaping is performed by operational 
amplifiers using input and feedback networks consisting of resistors and 
capacitors. Inherently, such analog circuitry is subject to drift, noise and 
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component aging. Replacement of the analog servo amplifier with digital tech- 
niques will improve reliability and reduce maintenance requirements. 

The block diagram of the analog servo amplifier using Laplace notation 
is shown in Figure 3.1. In order to implement these characteristics into a digital 
computer program it is necessary to convert the Laplace expressions into the 
Z -domain. To ease computational complexity, the block diagram is redrawn in 
Figure 3.2 with a slightly different arrangement leaving dynamic characteristics 
unchanged. Next a sampling frequency has to be established which is limited by 
the computer speed on one extreme and the computational accuracy on the other. 
The tracking loop has a bandwidth of approximately 1.0 Hz, however substantial 
attenuation is realized only after the first natural frequency of the structure 
which is around 5 Hz. Hence a sampling frequency of 20 Hz has been chosen. 

Conversion to the Z-domain follows a procedure outlined in reference 3.2 
(page 465 and following). This procedure makes use of the W-transform where 


w = 


z - 1 
Z + 1 


v 


t an 


ojT 


where 

a) = any frequency in S-domain, rad/sec 
v = corresponding frequency in W-domain rad/sec 
T = sampling period, sec 

The W-transform of A(s) can be found by making use of equation 3.2: 


(3.1) 

(3.2) 


A (s) 


1.96 

0.47 S + 1 



0.0532 rad/s 
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POSITION FEEDBACK 
(ENCODER OR SYNCHO) 


Figure 3-1. Block Diagram of Analog Servo Amplifier 


POSITION LOOP 



Figure 3-2. Modified Block Diagram of Analog Servo Amplifier 
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A (W) 


1.96 _ 1.96 

W l 18.9 W + 1 
0.0532 + 


Substitution of equation (3.1) into above expression yields the desired Z-transform 


A (Z) = 


1.96 

7-1 

18.9 - — - + 1 
Z + 1 


0.099 Z + 0.099 
Z - 0.899 


(3.3) 


In similar fashion B(s) and C(s) can be converted: 


B (s) = 0.75 


0.0864 s + 1 
0.47 s + 1 


and 


B (Z) = 


0.165 Z - 0.089 
Z - 0.899 


(3.4) 


C (s) = 12.1 


1.1 s + 1 

s 


C (Z) 


13.62 Z - 13.0 
Z - 1 


(3.5) 


Equations (3.3), (3.4) and (3.5) appear in the general form of: 


which lends itself readily for computer implementation. If a signal E. ^ (nT) 
is applied to G(Z), then the output E out (nT) is given by: 


E out ( n T ) = a o f ( n T) + f [(n - 1) T] 
f (n T) = E . (n T) - b. f [(n - 1) T] 

in 1 
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where 


T = sampling period, seconds 


n = number of sampling interval 
= 0, 1, 2, 3 . . . 

Figure 3.3a shows the step function response of the antenna using the analog 
servo amplifier and digital compensation. Figure 3.3b shows stepfunction 
responses using digital compensation for various sampling rates. Note that 
the selection of 20 samples per second is satisfactory. 

The digital compensation scheme is based on investigations reported in 
reference 3.3. 


3.2 Antenna Position Program 


In conventional tracking stations the Antenna Position Programmer (APP) 
is used to point the antenna in accordance with satellite predict data. This is 
accomplished by accepting predicts from the ENV to XY conversion program 
(covered in paragraph 3.3) interpolating the data and subtracting the antenna 
encoder reading. The resultant signal is an error which is fed to the servo 
amplifier. Whenever the antenna operates in this configuration, it is said to 
be in program mode. The APP is really a special purpose digital machine that 
accomplishes the tasks outlined. 

The function of the APP can be replaced by a program under the computer 
controlled antenna concept. Figure 3.4 shows the basic block diagram of the 
position programming function. The program accepts predict data at a one per 
second rate and synchronizes it with real time. The operator biases real time 
from the control console. Synchronized data is linearly interpolated to 20 data 
points per second. The operator can add position bias to this data. Subsequently 
the data is checked against physical limits of the antenna mount. Then the 
difference is formed between this quantity and the encoder reading. The resultant 
error is checked for its maximum limits and fed to the servo amplifier (or 
digital compensation program). 

Mathematically these relationships can be expressed as following: Time 
Bias: 




(n + T\ • V 

v bias-' 


n = 0, 1, 2,. • • 
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to 



Figure 3-4. Block Diagram of Antenna Position Programming 














where 


X 

X 


pn 

1 

pn 


^ias 

Interpolation: 


= synchronized predict value at time interval n 
= synchronized predict value modified by time bias 
= time bias introduced by operator 


X = X 1 

pm pn 


/ X p(n+1) " X pn\ 

\ 20 ) m T ’ m = °* l> 
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where 


T = sampling interval = 1/20 second 


X = predict value at time n + ^ 

pn> 20 


Position Bias: 


where 


x p m - ^ + x bii 


Xjj j as = position bias introduced by operator 


X om = predict value at time n + modified by position bias 


pm 


Servo Error: 


where 


E xm - X pm “ X em + X corr 


x D 


encoder reading at time interval n + 


m 

20 


X C orr = constant for encoder correction 


E xm = servo error at time interval n + ^ 

xm 20 
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In addition to the task described above this program drives the console 
displays as indicated in Figure 3.4. A history record is also generated which 
permits recording of such quantities: servo error, antenna position and velocity, 
age values, and time. The history record can be used for post-pass analysis. 


3.3 ENV to XY Conversion 

To aid satellite acquisition and tracking, prediction information has to be 
made available to the antenna system. Orbit calculations for satellites are 
made on the centralized computer complex at GSFC. Prediction information is 
transmitted from here to the various tracking stations. However, to save on 
transmission volume, local processing of prediction data is necessary which 
involves coordinate conversion and interpolation. This manipulating of data is 
handled by the PB250 computer at tracking station (UNIVAC 1248 for USB 
antennas). Under the Computer Controlled Antenna concept this function is 
handled off-line by the antenna system computer. 

The position of a satellite is originally determined using orbital elements 
in an E‘N ! V S (East, North, Vertical) coordinate system. The V'-axis is perpen- 
dicular to the mean earth equatorial plane and points to the north celestial pole. 
The E ! and N ! axes lie in the mean earth equatorial plane. The N T axis points 
toward the mean vernal equinox and the E 3 axis is situated to complete a right- 
hand coordinate system. The origin of this system is at the center of the earth. 
This system is then rotated and transformed into a topocentric ENV coordinate 
system with respect to the particular tracking site. In this system the V axis 
is directed along the upward vertical; the E axis along the eastward horizontal; 
and the N axis along the northward horizontal. 

The ENV vectors for a specific satellite are annotated with time and 
transmitted to the tracking station. Transmission is limited to a maximum of 
75 coordinate values for a particular pass of a particular satellite. This limi- 
tation also establishes the spacing of coordinate values in time. This spacing 
can vary from 1 to 20 minutes depending on the satellite orbit. Figure 3.5 shows 
a typical record that has been transmitted to a station. Predicts are generated 
for two weeks in advance and are transmitted to the station. 

In the Computer Controlled Antenna concept predicts as shown in 
Figure 3.5 are accepted and converted to predict records compatible with the 
antenna position program. The first step in this process is the conversion of 
ENV data to XY angular data which is accomplished using the following 
relationships: 
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TWX Address'! 
For NTTF j' 

Julian Day- 
Satellite ID 1 
Number (Nimbus 3)J 

Station Code] 
For NTTF J 

First & Last Date 
For Predict Period V 
(Year, Julian Day)J 

Date of Particular 
Pass (Year, Jul. Day) 



Number of Predict 
Points 


GTFT GWWW 
222 700810 
650370 1 
NfcTEST 
70208 2355 
70223 2255 




70222 15 5 ICC 


Special Data Terminal (Greenbelt, Md.) 

,Date (Year, Month, Day) 

Time of Pass (Hour, Minutes, Seconds) 

of Start & Finish of/ 

Period (Hour Minutes) 

Start of Pass 

(Hours, Minutes, Seconds) 
^Spacing Between Data 
Points (Minutes) 



Time 

Predict 


oc i-« 

00 10562-00 4 e 155-0007287 
C009880— 004451 1—0004457 
000877 1-004 l 45 5— CCC 1681 
0007840-0037534 0CC0524 


2 7 C 
275 
262 
254 


00054 B9-TT03 A 238 

"00X76 84 

”2 47 

0G05323— 002042C 

00C461 2 

22£ 

0004 146 — 0 0 2 6 452 

0 C C62G2 

225 

0002968-0022468 

CCC7746 

261 

0001793-0018361 

0CC8S40 

2 6 C 

0000627-0014186 

ccosee l 

2 6 1 

-00005 22-080593 7 

TCio^er 

274 

-ooaa 651-0005687 

00 10586 

272 

-0002753-000 1353 

00 i 1 148 

246 

— C00382 2 000251 1 

CC 1 l 048 

2 2 2 

-0004853 0007212 

CC 10666 

242 

— C005840 0011453 

CC 10C64 

246 

-0006778 0015740 

CC09 184 

26 7 

-0007664 0019535 

CCC6C5C 

267 

-0008492 0024073 

0CC6664 

2 5 1 

-0009259 0028130 

C C CS C J 2 

250 

-0009959 0032C54 

0CC3 162 

252 

-0010591 0085552 

CCC 1 0b8 

2 4 <4 

-001 1 183 0U3988C- 

-OOC7 271 

249 

— 00 11646 0043257- 

-CCC J820 

2 5 6 

-0012068 0046755- 

- CCC6S8 1 

268 


-25/ C648C> 


Z 

tNO 


PASS 


Data Points 

1 . Column: East Vector 

2. Column: North Vector 

3. Column: Vertical Vector 

4. Column: Checksum 

(All Vectors in KM Times 10) 


'Orbit Number 


Figure 3-5- Typical Predict Data Format as Sent to a Tracking Station 



£ 

X = arc tan — 
V 


Y = arc tan — ______ 

/E 2 + V 2 

Some additional conversions are performed which are useful in describing 
the particular pass to the operator. These are range, R; range rate, R; azimuth, 
A; and elevation, E. 


R = /E 2 + N 2 + V 2 


R ~ 




e sen t 


Spr, 


Time interval 


A = arc tan — 

N 

E= arc tan— — — 

/e 2 + N 2 

USB antennas are provided with an acquisition message that provide values 
for x and y angles directly, hence above conversions are not necessary. Using 
computed or provided X and Y angles with a period, T Q , between subsequent 
values, the data is interpolated to yield a fixed predict rate of one per second. 
The method employed is based on the Lagrange Five Point Interpolation Formula 
(Reference 3.4; Formula 25,2,15) which is given by: 


X t = ^ ( t2 - 1) t (t - 2) a 0 - i (t - 1) t (t 2 - 4) aj 
+ \ ( t2 - 1) (t 2 - 4) a 2 - i (t + 1) t (t 2 - 4) a 3 


+ 24 (t - 1} t (t + 2) a, 


(3.6) 
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where 


a 0 , a . 1 , a 2 , a 3 , a 4 = interpolation coefficients 
t = normalized time 
x t = angular value at time t 

Interpolation about a specific predict point, X q = 0 , is accomplished by evaluating 
the interpolation coefficients first using X q = 0 and two points on either side of 
X q = o. Evaluation of equation (3.6) at t = -2, -1, 0, 1, 2 and equating appropriate 
terms to predict values 


X t=-2 a 0 X q=-2 


X t=-1 = a i = V> 

X t=o a 2 X q=o 


X„_„ = = X „ 

t= 2 4 q= 2 


(3.7) 


Rearranging equation (3.6) gives: 


where 


x t = b 0 + b l 1 + b 2 t2 + b 3 t3 + b 4 t4 


b o = 


(3.8) 


b i = 


12 3 ° 3 3l + 3 3 ’ ' 


12 


i 1 2 5 2 1 

b = - — - a + _ a a. + — a. a A 

2 24 0 3 1 4 2 3 3 24 4 


b 3 = 


12 


a + _ a, a, + — a A 

0 6 1 6 3 12 4 


b 4 = 


11111 

24 3 ° 6 3l + 4 32 " 6 83 + 24 84 
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Coefficients b 0 through b 4 can be numerically evaluated using relationships in 
(3.7). Predict angles at a one per second rate can now be generated by subsequent 
solutions of equation (3.8) for different values of time, t, such that: 


t = 


n 1 _ 

60 T "T 
q 


n = 0, 1, 2. . . (60 T q - 1) 


where T q = spacing between predict values in whole minutes (1 < T q ^ 20). 

After computation of 60 T q predict values new interpolation coefficients 
are established, by dropping the last predict value, X q = _ 2 , and using new 
point X q = 3 . This process is then continued until the whole pass has been 
covered. 


3.4 Scan Pattern 

To aid in the acquisition phase antennas are provided with a scan generator. 
This generator superimposes a search pattern in most modes of operation. The 
pattern used most frequently is that of an Archimedes spiral. Conventionally, 
implementation of this pattern is accomplished by using electromechanical devices 
such as a sine-cosine pot or a resolver and an electrical or mechanical integrator. 
In the computer controlled antenna the equation for a spiral is solved and then 
added to the current X and Y angles of the antenna. The equation for a spiral 
is given by: 


X = a t cos a> t 

Y = a t sin w t (3.9) 


where 


a = spiral width 

oj = spiral angular frequency. 

The difficulty with this equation is that the tangential velocity changes. 
This change is undesirable when one investigates acquisition probabilities but 
is difficult to overcome with the electromechanical scan generator. In the 
computer controlled antenna concept this shortcoming can be eliminated. The 
tangential velocity, v , can be kept essentially constant by making a substitution 
into equation (3.9): 
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X - a z cos co z 


Y = a z sin co z 


(3.10) 


where 


The tangential velocity is given by: 


Z = f (t) 


v = /x 2 + y 2 


(3.11) 


The derivatives of equation (3.10) are: 


X = a z co 


s oi z r a at z z sinojz 


Y = a z sin 


w z + a w z z cos co z 


substitution of these derivatives into (3.11) yields: 


= a z / 1 + co 2 z 2 


and 


for values of z such that 


v a a o) z z 




2 z 2 > > 1. 


A function that satisfies v = constant is 


z = /t and z = — — L. 

2 /t 


If the spiral is adjusted to 2-1/2 turns then 


co Z = 5 77 
max 
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and defining a maximum radius, r max : 


r - az 

max max 


it is possible to find expressions for constants a and co in equation (3.10): 


= 0. 36 / v r 


co = 5.6 


Hence the scan pattern implemented on the computer is: 


X = O' 36 /v r max 11 COS 


5.6 


/V t 


Y = 0.36 /v r max t sin 


5.6 




The operator specifies, r , by setting the scan radius switch on the operators 
console, while the tangential velocity, v, is fixed at one degree per second. It 
will be substantially simple to extend this algorithm to one that maximizes 
acquisition probability by automatically selecting v and r max based on radio- 
frequency receiver selection. Figure 3.6 shows the scan pattern generated by 
the above method. 


4. AUGMENTATION 

The augmentation mode is a combination between autotrack and program 
mode and displays a tracking error that is less than either of the two other 
modes. In the augmentation mode use is made of the inherent characteristics 
of the autotrack and program modes: 

Autotrack: Major tracking error due to receiver noise 

Program: Major tracking error due to timing inaccuracies and misalign- 

ments of antenna radio frequency axis. 
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In order to analyze this mode of operation a satellite in a 160.9 km (100 mile) 
circular orbit is assumed passing directly over the tracking antenna. Tracking 
errors in the various modes of operation will be evaluated and compared with 
each other. 


4.1 Satellite Motion 

In order to analyze the tracking performance of the antenna equivalent sinusoidal 
inputs will be used. These inputs are computed using a technique described in 
reference 4.1 (page 43 and following). Use of sinusoidal inputs simplifies compu- 
tational manipulations. The orbital velocity of the satellite is established using: 


V 



where 


v = velocity of satellite, km/s 
k = gravitational constant 
= 398,297.83 km 3 /s 2 
r = semimajor axis 
= earth radius + satellite height 
= 6377.8 + 160.9 = 6538.7 km 


Hence 


7398.297 
V 6538.7 


7.805 km/s 


A normalization constant, a, can be computed 


_ v _ 7.805 
h 160.9 


0.0485 1/s 


where 


h = satellite height = 160.9 km 
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The equivalent sinusoid for antenna position is given by: 

9 (t) ; A sinw t 

Sp v 7 p Sp 

where 

<9 sp (t) = amplitude of equivalent sinusoid due to position changes, rad 

A = peak amplitude of equivalent sinusoid due to position changes, 
P rad 

o) sp = frequency of equivalent sinusoid due to position changes, 
rad/s 

From a curve given in reference 4.1 (page 46) 

A p = 0.64 x 2 = 1.28 rad 

oj sp = ll a x 0.0485 = 0.0254 rad/s 

Similarly, using subscripts v, a, and j for velocity, acceleration, and jerk 
respectively 

A v = 0.5 a = 0.0243 rad/s 

co = a = 0.0846 rad/s 
sv 3.6 

A a = 0.65 a 2 = 0.00153 rad/s 2 

co = ^ 77 3 - = 0.152 rad/s 
s 2.0 

A. = 0.6 a 3 = 0.0000685 rad/s 3 
"sj = ~Fo~ = 0-305 rad/s 
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Hence the equivalent sinusoids for the stipulated satellite are: 


0 (t) = 1.28 sin 0.0254 t 

6 sv (t) = 0.0243 sin 0.0846 t 
9 sa (t) = 0.00153 sin 0.152 t 
0 sj (t) = 0.0000685 sin 0.305 t 


4.2 Servo Loop 


Before discussing a particular mode of operation it is necessary to analyze 
the servo loop without considering receiver, prediction and alignment errors. 
As was shown in paragraph 3.1 the basic servo loop is a type II servo. Such a 
servo can be stabilized by a simple lead-lag network and has an open loop 
transfer function G : 


K a (T, s + 1) 
G = 

SO 

s 1 2 (T 2 s + 1) 

where 


K a = acceleration gain, 1 /seconds 2 
T x = time constant of lead, seconds 
T 2 = time constant of lag, seconds 

For a phase margin of 0.662 rad and a damping factor of 0.5 the following re- 
lationships can be established with the crossover frequency, u> : 

C 


K 


a 


1 

2 


O) 


2 

c 


1 

T ! 



2 



2 w 


C 
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Hence the open loop transfer function can be rewritten: 


- (-L s + 1 

g = : 

30 2/1 A 

S ' 4 f s + 1 

\ 2 " c > 


The closed loop transfer function, G , can be found: 


s + 1 


G = — 

s 1 + G. 


’ 1 A / 1 2 1 

s + 1] / s 2 + — s + 1 


Similarly the error transfer function can be established: 


, \ s h— s + M 

6 S ^ _ 1 _ 2 \ 2 "c / 

<*> ' iTgt. = s j + J. s + 1 


where 


e s = servo error, rad 
0 - command angle input, rad 


Assuming sinusoidal excitation: 


6 ( t) = A s in a) t 


S (s) = 


u - 

s 1+J_S 2 


the error in equation (4.3) can be expressed as: 
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c / \ 2 A 

e s (S) = 

0) 2 "s 
c 


s 2 


2 gj 


s + 1 


1 s 2 +_i s + i\ 41 s + l 


. co* 


00 




1 .2 


s 2 + 1 


oo A 


Using partial fraction expansion e (s) can be written as: 


e (s) = co A 


‘l 


l 


CO CO CO CO 

S + — + j — /3 s + — - j — /3 
2 2 2 2 


k 2 k 3 k 3 

+ + + 


s+w s+ioj s - ] u 

c J s J s 


where the * indicates the complex conjugate. 

By inspection, only k 3 terms will yield a steady state response, hence 


e (s) I 

8 steady state 


S + J Cd s - j Cd 


sj 


co, A 


and 


k 3 = 


S 2 (S + 2 co ) 


(s 2 + co s + co 2 ) (s + co ) (s + j co ) 


s=j to- 


co 


('P 


2\l/2 


j'P i 


L" 2 (CO 2 - 2 a ?) 2 + CO 2 (co 2 ~ 2 W 2 )] 1/? e 3 ^ 


where 


(4.4) 
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tan 



tan 


( M c ~ 2 
("s ~ 2 "c) 


0 = 0i - 02 

Transformation of equation (4.4) into the time domain yields: 


€ (t) 

s steady state 


(4 


d ) 2 + £0 2 ) 1/2 

C S'' 


W - 2 co 2 ) 2 + col (col - 2 0)2)2] 1/2 


cos 


(a) t + </<) 


and for the case co <<co : 

S C 


€ 

s 


(t)| 


steady state 



COS (o) g t + l/>) 


(4.5) 


Amplitude A and frequency co s of the excitation have been established in 
paragraph 4.1. There remains the crossover frequency, o> c , to be found. This 
can be done following the procedure in reference 4.1. An approximate transfer 
function established by satellite motion can be found using the results from 
paragraph 4.1 and is shown graphically in Figure 4.1. Adjusting the gain of 
transfer function G such that the approximate transfer function is below 

SO. 

G so and = 6.28 rad/s is found. This indeed represents the design value of 
most STDN antennas. It should also be noted that the design as chosen provides 
sufficient attenuation for the structural natural frequency of the antenna which 
occurs around 5.0 Hz. / 

The maximum steady state errors can now be computed using equation (4.5) 


0 0 2 54 2 

e /position = 2 — (1.28) = 0.0419 m rad 

6. 28 2 


e /velocity = 2 0 08462 (0.0243) = 0.00879 m rad 
S 6. 28 2 
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. (ERROR = 0.055m rad) _ 
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\ 

JERK © 
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Figure 4-1. Bode Plot for Tracking Loop 


V 




e g /accel er at ion = 2 


(0.00152) = 0.00178 m rad 

6.28 2 


0 305 2 

e s /j erk = 2 -11)1... (0.0000685) = 0.000323 m rad 
6. 28 2 

As can be seen from above the major servo error is due to the equivalent sinusoid 
of position. 


4.3 Autotrack Mode 

Next the operation of the antenna in the autotrack mode will be considered. 
Operation in this mode is essentially identical to that of the closed loop servo 
analyzed in the previous paragraph. The only difference is a receiver noise 
term that has to be added to the loop. But before discussing the autotrack loop, 
the pertinent receiver characteristic will be established. The assumption is 
made that the noise power spectrum, <P TT , is flat up to a frequency, , which is 
the bandwidth of the receiver error channel and that the spectrum has an 
amplitude K R . Hence the spectrum can be expressed mathematically: 


= K * 


to 2 + CO 2 


(4.6) 


If the additional assumption is made that the root-mean-square noise of the 
receiver, E (rms), has a value of 2 mrad and that the error channel bandwidth 
is 62.8 rad/s (ten times servo bandwidth) then K R can be evaluated: 


and 


RR 


( rms) 


-J 

Jo 


^RR O) d W = 


77 k R W R 




2 S RR ( rmS V 


n 


2 x 4 = 0.00404 m r ad 3 / S 

62.8 77 
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The tracking error in the autotrack mode can be found by applying the 
superposition theorem. The tracking error is the sum of the error due to 
satellite motion alone and the error due to receiver noise alone. The first 
error is the servo error, E s , established in the previous paragraphs while the 
second error component, Eg , will be developed below. Using reference 4.2 as a 
basis: 


e 


2 

R 


1 

2 tt j 



l G s (s)| 2 (s) d s 


G s (s) is the closed loop transfer function developed previously and is given in 
equation (4.2). The spectrum $ (s) is the transform of spectrum c/> RR (w) given 

by (4.6). Hence above expression can be rewritten as: 


S 2 = k R -c 4 -4 


1 1 

- j® 

2 s co 

C 

2 

d s 

tO! 

3 

-j 00 

s 3 + 2 co s 2 + a; 3 

C c 

2 2 
5 +a) l 


This expression can be put into a form compatible with reference (4.2) (Appendix E): 


e R = k R < 




c (s) C (- S) 
d (s) d (- s) 


d s 


(4.7) 


where 

c (s) = 2 s + co c 

d (s) = s 4 + (2 w c + cj r ) s 3 + 2 «. (co c + u) R ) s 2 
+ OJ 2 (co e + 2 CO R ) S + OJ 3 c OJ R 
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Making use of a previous assumption that co = 10 co the coefficients of the 

K c 

above expression are: 

C 0 = "c d 0 = 10 "c 

Cj = 2 dj = 2 1 

d 2 = 22 co\ 

d 3 = 12 

d 4= 1 

From reference (4.2): 

j = C l d 0 d 3 d 4 + C 0 (- d l d 4 + d 2 d 3 d 4> 

2 d 0 d 4 d 0 d 3 " d l d 4 + d l d 2 d 3> 

Substitution of the coefficients and simplification yields 

I 4 = 0.00988 oT 5 

and substitution of this quantity into (4.7) gives 

4 = °'988 Kjj 

and 

e R = 0.995 /k R w c = 0.99 5 /(0.00404) (6.28) 

= 0. 158 mrad 


4.4 Program Mode 

The major errors in the program mode of operation are due to prediction 
errors and misalignments fo the radio frequency axis of the antenna. Both of 
these errors are of the dc bias type and therefore are not affected by the opera- 
tion of the servo. To express this phenomenon mathematically, it is only 
necessary to consider an input to the servo loop that is the sinusoid due to 
satellite motion, but shifted by a bias term which is proportional to the error. 
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The transfer function of the programmer which includes antenna misalign- 
ments can be computed by noting that the input is 

0 s (t) = Asinw t 


and the output is 


9 p (t) = A sin (co s t + <p) 

where 

cp = phase shift proportional to error term. 

Since 4> is generally small, sin 4>= 4> and cos <p « 1. With this 8 (t) can be 

expressed: 


9 p (t) = A (sin w s t + cfr cos a> s t) 

Taking Laplace transforms: 


0, (s) 


A 


1 + 


CO s 



s 


4> 


s + 1 


1 + — - — s 2 
co s 2 _1 


Hence the transfer function, G p (s), can be found: 


G p( s)=^lL± stl 

P 8 (s) CO 

S'' R 


(4.8) 


If the assumption is made that program error, e p , is equal to autotrack error, 
then 


6 P = 9 S (*) - e t (t) = - 4> A COS o) s t 

e n (t)| • = 0.158 m r ad 

P max 
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The phase shift, 0 , can be numerically evaluated for the equivalent sinusoid of 
position: 


0 = x io -3 = 0.000158 rad 

1.28 


4.5 Augmentation Mode 


The basic approach to this mode of operation is based on an idea and pre- 
liminary investigation by G. C. Winston (reference 1.1). In this mode the antenna 
is operated as if it were in the program mode, but the predicted angle is modified 
by a term derived from the receiver. This correction term is obtained by passing 
the receiver error through a low pass filter and then performing an integration. 
Figure 4.2 shows the basic block diagram of this mode. The augmentation filter 
has a transfer function: 


G f (s) 


s (T f s + 1) 


where 


k f = gain of filter 
T f = time constant of filter, s. 
The constants will be adjusted such that 


Hence, 


10 K, = co and 

t C 


T f = — 

f Kx 


Gf (s) = 


CO 

c 

To~ 




f 10 


s + 1 


s (s + K f ) 




(4.9) 


Analysis of this mode will make use again of the superposition theorem, i.e., 
first errors will be computed with receiver noise input only (# s (t) = 0) and then 
with satellite motion only (<£rr(") = 0). The effective transfer function, G R , to 
which receiver noise is applied is: 
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Figure 4-2. Block Diagram for Augmentation Mode 


G f (s) G (s) 

n (s) = 

^ 1 +G f (s) G s (s) 

From the development of expression (4.8) it can be seen that the bandwidth of the 
augmentation filter is approximately one tenth of the servo loop. Hence G R (s) 
can be simplified and expression (4.9) can be substituted: 


GrC 5 )- 


G { (s) 

1 + G f (s) 


Kf 

+ K f s + Kf 


Using, again, techniques from reference (4.2) the error due to receiver noise, 
e R , will be computed: 
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.12 __L_ r c 

R J-jc 


I G A C S )I 2 ^RR ( S ) d S 




1 

r }<» 

1 

2 ds 

2 77 j 

j CO 

s 8 + K f s + K* 

s 2 + "r 


Evaluating the integral and noting: 

&J R = 10 co Q - 100 k f 
the following expression is obtained: 


£* = 0.707 /Kjj K f 

Using numerical values for K and k : 

R i 

e x R = 0.707 /(0. 00404) (0.628) = 0.0356 mrad 

Next the error due to satellite motion will be investigated. Following the 
block diagram it can be shown that the transfer function 0 A /q s is: 

e K ( s ) _ G s O) G p (s) G s (s) Gp (s) 

& s CO ~ 1 + G s (s) Gp (s) + 1 + G s (s) Gp (s) 

Using the previous assumption this transfer function can be approximated as 


e K (0_ G, (s) Gp (s) G p (s) + Gp (s) 

e s (s) " 1 + Gp (s) + 1 + Gp (s) + 1 + Gp (s) 
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The error, e 1 , under these conditions is: 


e p ( s ) = 6 s ( s ) - 9 k (S) 


1 " G P „ 

9 (s) 

1 -Gp 00 s 


and substitution of equations (4.8) and (4.9) yields: 


0 „2 


CO 


(s + K f ) 


«;c s)=-f r (») 

s 2 + s Kf + Kf 


Applying sinusoidal excitation and using partial fraction expansion: 


■^ s2(s + K f ) A a) 

1 . v s s 

( S ) = 

P v v o 2 2 

-f-2 + oo* 


S 2 ' + S K f + K 2 s 2 


K, 


K 


K f 


v f . /3 _ K c . /3~„ 

s+ y +J — K f S+ T- J ~ 


K„ 


s + J " s S - J W s 


By inspection only K 2 terms will yield a steady state response hence: 


e 1 (s)|- Ka ! K 2 

P steady state s + i CO S - i 

J S J t 


where 
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E p (t)| = 0.00632 mrad 

** max 

With this information it is possible to compare errors in the various modes: 
Autotrack: e auto = ^ e s + £ r " mrad 

Program: e prog = ^ + e l ~ 0163 mrad 

Augmentation: e aug = l/e s + £ r 2 + £ p 2 = 0.0553 mrad 

as is evident, tracking in the augmentation mode improves the accuracy by a 
factor of three. 

It should be recalled, however, that the assumption has been made that 
program error and autotrack error are equal. For situations in which the 
program error is larger, due to severe misalignment of antenna boresite for 
example, the factor of improvement of the augmentation mode increases considerably. 
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Figure 4.3 shows the performance of the augmentation mode while tracking 
Nimbus IV at 136 MHz. Initially the antenna operates in the autotrack mode. 
Then the operator switches to program mode. Note the off-set in this mode 
which is due to antenna misalignments and predict errors. Finally the augmen- 
tation mode is selected and the error goes to essentially zero after one overshoot. 
While in program mode the operator could have introduced angular bias to reduce 
the off-set error, however this bias varies as the satellite pass progresses and 
is automatically achieved by selecting the augmentation mode. 
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Figure 4-3. Performance of the Antuentation Mode 


5. SUMMARY 

The computer controlled antenna system has several advantages over the 
conventional system: reliability and consistency of operation, minimization of 
turn-around time between satellite passes, reduction on floor space require- 
ments, fast identification of component or subsystem failures or misadjustments 
and in the case of new antenna procurements substantial savings on special 
equipment requirements. But above all there is one additional attribute and 
that is flexibility. Only software changes are required to change tracking 
algorithms, data to be monitored and to change the degree of automation desired. 

All software programs have been written and checked using the 40 ft. antenna 
and the XDS Sigma 5 computer at NTTF. Test results have been extremely 
encouraging. Currently, ah effort is in progress to duplicate this system using 
a smaller dedicated computer, a DEC PDP-ll/20, in conjunction with the 40-ft. 
antenna. This new system will represent the final prototype for STDN antennas. 
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